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DETECTING BINOMIALITY 


CARSTEN CONRADI AND THOMAS KAHLE 


Abstract. Binomial ideals are special polynomial ideals with many algorithmically 
and theoretically nice properties. We discuss the problem of deciding if a given poly¬ 
nomial ideal is binomial. While the methods are general, our main motivation and 
source of examples is the simplification of steady state equations of chemical reaction 
networks. For homogeneous ideals we give an efficient, Grobner-free algorithm for 
binomiality detection, based on linear algebra only. On inhomogeneous input the algo¬ 
rithm can only give a sufficient condition for binomiality. As a remedy we construct a 
heuristic toolbox that can lead to simplifications even if the given ideal is not binomial. 
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1. Introduction 

Non-linear algebra is a mainstay in modern applied mathematics and across the 
sciences. Very often non-linearity comes in the form of polynomial equations which 
are much more flexible than linear equations in modeling complex phenomena. The 
price to be paid is that their mathematical theory—commutative algebra and algebraic 
geometry—is much more involved than linear algebra. Fortunately, polynomial systems 
in applications often have special structures. In this paper we focus on sparsity, that is, 
polynomials having few terms. 

The sparsest polynomials are monomials. Systems of monomial equations are a big 
topic in algebraic combinatorics, but in the view of modeling they are not much help. 
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Their solution sets are unions of coordinate hyperplanes. The next and more interesting 
class are binomial systems in which each polynomial is allowed to have two terms. 
Binomials are flexible enough to model many interesting phenomena, but sparse enough 
to allow a specialized theory [8]. The strongest classical results about binomial systems 
require one to seek solutions in algebraically closed held such as the complex numbers. 
However, for the objects in applications (think of concentrations or probabilities) this 
assumption is prohibitive. One often works with non-negative real numbers and this 
leads to the helds of real and semi-algebraic geometry. New theory in combinatorial 
commutative algebra shows that for binomial equations held assumptions can be skirted 
and that the dependence of binomial systems on their coefficients is quite weak [14]. For 
binomial equations one can hope for results that do not depend on the explicit values 
of the parameters and are thus robust in the presence of uncertainty. 

The main theme of this paper is how to detect binomiality, that is, how to decide if a 
given polynomial system is equivalent to a binomial system. The common way to decide 
binomiality is to compute a Grobner basis since an ideal can be generated by binomials 
if and only if any reduced Grobner basis is binomial [8, Corollary 1.2]. For polynomial 
systems arising in applications, however, computing a reduced Grobner basis is often too 
demanding: as parameter values are unknown, computations have to be performed over 
the held of rational functions in the parameters. Even though this is computationally 
feasible, it is time consuming and usually yields an output that is hard to digest for 
humans. This added complexity comes from the fact that Grobner bases contain a 
lot more information than what may be needed for a specihc task such as deciding 
binomiality of a polynomial system. Hence Grobner-free methods are desirable. 


Grobner-free methods. Grobner bases started as a generalization of Gauss elimina¬ 
tion to polynomials. They have since come back to their roots in linear algebra by 
the advent of F4 and F5 type algorithms which try to arrange computations so that 
sparse linear algebra can exploited [7]. Our method draws on linear algebra in bases of 
monomials too, and is inspired by these developments in computer algebra. 

Deciding if a set of polynomials can be brought into binomial form using linear alge¬ 
bra is the question whether the coefficient matrix has a partitioning kernel basis (Dehni- 
tion 2.1 and Proposition 2.5). Deciding this property requires only row reductions and 
hence is computationally cheap compared to Grobner bases. It was shown in [18] that, 
if the coefficient matrix of a suitably extended polynomial system admits a partitioning 
kernel basis, then the polynomial system is generated by binomials. As a hrst insight 
we show that the converse of this need not hold (Example 2.8). 

In general computer algebra prohts from homogeneity. This is true for Grobner bases 
where, for example, Hilbert function driven algorithms can be used to convert a basis 
for a term-order that is cheap to compute into one for an expensive order, such as lex. 
We also observe this phenomenon in our Grobner-free approach: a satisfying answer 
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to the binomial detection problem can be found if the given system of polynomials is 
homogeneous. In Section 3 we discuss this case which eventually leads to Algorithm 3.3. 

In the inhomogeneous case things are more complicated. Grobner basis computations 
can be reduced to the homogeneous case by an easy trick. Detection of binomiality can 
not (Example 4.1). We address this problem by collecting heuristic approaches that, in 
the best case, establish binomiality without Grobner bases (Recipe 4.5). Our approaches 
can also be used if the system is not entirely binomial, but has some binomials. In 
Example 4.4 we demonstrate this on a polynomial system from [5]. 


Binomial steady state ideals. While binomiality detection can be applied to any 
polynomial system, our motivation comes from chemical reaction network theory where 
ordinary differential equations with polynomial right-hand sides are used to model dy¬ 
namic processes in systems biology. The mathematics of these systems is extremely 
challenging, in particular since realistic models are huge and involve uncertain param¬ 
eters. As a consequence of the latter, studying dynamical systems arising in biological 
applications often amounts to studying parameterized families of polynomial ODEs. 
The hrst order of business (and concern of a large part of the work in the area) is to 
determine steady states which are thus the non-negative real zeros of families of pa¬ 
rameterized polynomial equations. Moreover, the structure of the polynomial ODEs 
entails the existence of affine linear subspaces that are invariant for solutions. Hence 
questions concerning existence and uniqueness of steady states or existence of multiple 
steady states are equivalent to questions regarding the intersection of the zero set of a 
parameterized family of polynomials with a family of affine linear subspaces. 

If the polynomial equations describing steady states are equivalent to binomial equa¬ 
tions (that is, generate a binomial ideal), then their mathematical analysis becomes 
much easier. This is the main theme of [18]. If a system is binomial, then, for in¬ 
stance, one can decide efficiently if positive steady states exist. If so, then a monomial 
parametrization can be found using only linear algebra over the integers [8, Section 2], 
and the steady states are toric. they are the positive real part of a toric variety. A 
sufficient criterion for toric steady states appears in [13, Theorem 4.1]. Since zero sets 
of general polynomial systems need not have parametrizations at all, we view the task 
of detecting binomiality as an important step in analyzing polynomials—in systems 
biology, or other areas like algebraic statistics, control theory, economics, etc. 

A frequent and challenging problem in the analysis of dynamical systems in biology 
is to decide multistationarity, that is, the existence of parameter values leading to more 
than one isolated steady state. A variety of results for precluding multistationarity has 
appeared in recent years. See, for instance, [12, 21, 3, 1] for methods employing the 
Jacobian. Similarly, several sufficient criteria for multistationarity have emerged (for 
example [4, 5]). In general this problem remains very hard. However, in the case of bi¬ 
nomial steady state equations, the question of multistationarity can often be answered 
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effectively, for example positively by [18, Theorem 5.5], or negatively by [17, Theo¬ 
rem 1.4]. Both of these results require only the study of systems of linear inequalities. 

Notation. In this paper we work with the polynomial ring ]k[a:i,..., x„] in n variables. 
The coefficient field k is usually R, or the field of real rational functions in a set of 
parameters. Our methods are agnostic towards the field. A system of polynomial 
equations /i = /2 = • • • = /^ = 0 in the variables xi,... ,Xn is encoded in the ideal 
{fi, ..., fs) C k[a:i,..., Xn]. A polynomial is homogeneous if all its terms have the 
same total degree, and an ideal is homogeneous if it can be generated by homogeneous 
polynomials. A binomial is a polynomial with at most two terms. In particular, a 
monomial is a binomial. It is important to distinguish between binomial ideals and 
binomial systems. A binomial system /i = • • • = /^ = 0 is a polynomial system such 
that each fi is a binomial. In contrast, an ideal (/i ,..., fg) is a binomial ideal if there 
exist binomials that generate the same ideal. Thus general non-binomials do not form 
a binomial system, even if they generate a binomial ideal. For the sake of brevity we 
will not give an introduction to commutative algebra here, but refer to standard text 
books like [6, 2]. The very modest amount of matroid theory necessary in Section 2 can 
be picked up from the hrst pages of [20]. 

Acknowledgment . We thank David Cox for discussions on the role of homogeneity in 
computer algebra. We are grateful to Alicia Dickenstein for pointing out a crucial error 
in an earlier version of this paper. TK is supported by CDS the Center for Dynamical 
Systems at Otto-von-Guerieke University. 

2. Grobner-eree criteria for binomiality 

The most basic criterion to decide if an ideal is binomial is to compute a Grobner 
basis. This works because the Buchberger algorithm is binomial-friendly: an S-pair 
of binomials is a binomial. Since the reduced Grobner basis is unique and must be 
computable from the binomial generators, it consists of binomials if and only if the 
ideal is binomial. However, Grobner bases can be very hard to compute, so other 
criteria using only linear algebra are also desirable. Linear algebra enters, when we 
write a polynomial system as A'^{x), the product of a coefficient matrix A with entries 
in k, and a vector of monomials ^/(x). Glearly, if we use row operations on the matrix 
to bring it into a form where each row has at most two non-zero entries, then the ideal is 
generated by binomials and monomials. This criterion is too naive to detect all binomial 
ideals since it allows only k-linear combinations of the given polynomials. We show in 
Section 3 that, at least for homogeneous ideals, it can be extended to a characterization. 
Before we embark into the details, we formalize the condition on the matrix. 

Definition 2.1. A matrix A has a partitioning kernel basis if its kernel admits a basis 
of vectors with disjoint supports, that is, if there exists a basis b^^\ ... of ker(A) 
such that supp(6^*)) fi supp(6(-^^) = 0 for all i ^ j. 
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The following proposition allows one to check for a partitioning kernel basis with 
linear algebra. The nnderlying reason is the very restricted strnctnre of the kernel, 
expressed best in matroid langnage. 

Proposition 2.2. The following are equivalent for any matrix A. 

1) A has a partitioning kernel basis. 

2) The column matroid of A is a direct sum of uniform matroids 17r-i,r of corank 
one, and possibly several coloops Ui^i. 

3) The reduced row echelon form of A has at most two non-zero entries in each 
row. 

Proof. 1 2: Let bi,... ,bk he the partitioning kernel basis. The snpports of this basis 

satisfy the circnit axioms and are thus equal to the circuits of the column matroid of A. 
Indeed, non-containment and circuit elimination are satished trivially because there is 
no overlap between any two circuits. For any non-zero element b G ker(yl), we have 
b = By the partitioning kernel basis property supp(6) = 7^ 0}) 

so either b is proportional to one of the bi, or its support properly contains the support 
of a circuit, so it cannot be a circuit. The columns of A which do not appear in any 
circuit are coloops and the remaining columns form a direct sum of k uniform matroids 
of corank one. 

2 3: If the column matroid of A is a direct sum of matroids, then the unique reduced 
row echelon form has block structure corresponding to the direct sum decomposition. 
Therefore it suffices to consider a single block which has one-dimensional kernel of full 
support (the coloops are (1 x l)-identity blocks). Ignoring zero rows, the reduced row 
echelon form of such a matrix is (1^-1 |c) where r — 1 is the rank and c G 

3 1: Rows of the reduced row echelon form with exactly one non-zero entry corre¬ 

spond to positions where every element of the kernel has a zero. Thus we can assume 
that there are none and each row of A has exactly two non-zero entries. Let c be a 
non-pivotal column with r — 1 non-zero entries. The restriction of A to c and the cor¬ 
responding pivotal columns yields a block containing {Ir-i\c) and some zero rows. The 
unique kernel vector corresponding to the dependencies in this block is orthogonal to 
the kernel of the remaining columns. This procedure can be applied to any non-pivotal 
column. The thus constructed basis is a partitioning kernel basis. □ 

Remark 2.3. Proposition 2.2 shows that the complexity of deciding if a matrix has 
a partitioning kernel basis is essentially the same as that of Gauss-Jordan elimination. 
One needs 0{n^) held operations where n is the larger of the dimensions of the matrix. 

Remark 2.4. A direct sum of (arbitrary) uniform matroids is called a partition matroid. 

We now translate Proposition 2.2 to polynomial systems. 

Proposition 2.5. If A has a partitioning kernel basis, and \k(a;) is a vector of mono¬ 
mials of appropriate length, then the ideal (AT(x)) C k[xi,...,x„] is binomial. If 
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(x) is any system that can be transformed into a binomial system using only k-linear 
combinations, then A has a partitioning kernel basis. 

Proof. Up to coloops the first part is [18, Theorem 3.3] and the coloops only give 
monomials. The second statement is clear since k-linear combinations of polynomials 
are row operations on the coefficient matrix and those do not change the kernel. □ 

Onr general strategy is to snitably extend a given system A^{x) with rednndant 
polynomials snch that Proposition 2.5 yields binomiality of an extended system A''^'{x). 
This happens in the following example. 

Example 2.6. Let fi = x — y, f 2 = z — w, and /a = x{fi + / 2 ) = x"^ — xy + xz — xw. 
Ordering the monomials = {x, y, z, w, x"^, xy, xz, xw), the system linearizes as 

/l-lO 0 0 0 0 0\ 

/ = (/i,/2,/3)= 0 0 1 -1 0 0 0 0 -T. 

\0 0 0 0 1 -1 1 -1/ 

The coefficient matrix is in rednced row echelon form and does not have a partitioning 
kernel basis by Proposition 2.2. Algorithm 3.3 takes this problem into acconnt, working 
degree by degree. The ideal is binomial since /a is a binomial (in fact, zero) in the 
qnotient ring k[x, y, z]/{x — y, z — w). 

The hrst hint into how to extend A and \k(a;) is the following theorem dne to Perez 
Millan et ah 

Theorem 2.7 ([18, Theorem 3.19]). Let fi = f 2 = ■ ■ ■ = fs = ^ be a polynomial system. 
If there exist monomials ...,a;"™- such that, for some ii,... ,im G [s], the system 

( 2 . 1 ) /^ = /2 = ... = /, = = ••• = = 0 

has a coefficient matrix with a partitioning kernel basis, then {fi,...,fs) is binomial. 

Theorem 2.7 is trne since the additional generators in (2.1) do not change the ideal 
that the system generates. This together with the explicit description of the binomial 
generators in the case of a partitioning kernel basis [18, Theorem 3.3] yields the resnlt. 
If the condition in Theorem 2.7 were also necessary, then a test for binomiality conld 
be bnilt on trying to systematically identify the monomials t". However, the converse 
of Theorem 2.7 is not trne. 

Example 2.8. Let I = (/i,/ 2 ) be the homogeneous binomial ideal generated by the 
non-binomials fi = x — y + x"^ + y‘^ + z‘^, f 2 = x‘^ + y"^ + z‘^. For no choice of monomials 
mil, ■ ■ ■ 5 f^ir, m 2 i,..., rn 2 s G k[a;, y, z] does the coefficient matrix of the system 

(2.2) fi = mii/i = • • • = mi,./i = /2 = ^ 21/2 = • • • = m2sf2 = 0 

have a partitioning kernel basis. 

For the proof Example 2.8 we hrst need the following curious little fact. 
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Lemma 2.9. The ideal I = {x^ + ?/^ + z^) C k[a;, y, z] does not contain a non-zero 
binomial. 

Proof. We can assume that k is algebraically closed, since if I contains a non-zero 
binomial, then so does its extension to the algebraic closure. Assume that for some 
/ G k[a:, ?/, z], the product b = f{x^ z"^) G I is a. binomial. We can assume that 

b is not divisible by any variable. Indeed, if a variable divides b, then it divides / and 
we find a lower degree binomial in I. Since I is homogeneous, we can also assume that 
b is homogeneous. Potentially renaming the variables, we can assume 6 = — Xy^z'^~^ 

for some 0 < s < d and A G k. Since k is algebraically closed, there is a solution ^ to 
the equation x^ = —1. The generator x^ -f--|- z"^ vanishes at (^, 1, 0) but b does not 
vanish there. This contradiction shows that I cannot contain a binomial. □ 

Proof of Example 2.8. Let d be the highest total degree among monomials in the sys¬ 
tem (2.2) and consider the restriction of all involved polynomials to degree d. Since the 
highest degree part of both /i and /2 equals (x^ -+■ y"^ -\- z^), only monomial multiples 
of (x^ -h y^ + z^) can contribute to the degree d part of the system (2.2). If the whole 
coefficient matrix of (2.2) had a partitioning kernel basis, then also the submatrix with 
only the columns for degree d monomials had one (Proposition 2.2). In this case row 
reductions on the submatrix would yield a binomial in degree d in the ideal generated 
by (x^ + y'^ + z‘^). This is impossible by Lemma 2.9. □ 

Example 2.8 may seem contrived, but this kind of “trivial obfuscation” of binomials 
does happen in applications. Of course, for humans it is obvious that one should first 
isolate the linear binomial x — y and then search for implied quadratic binomials which 
reduce the trinomial. Our next aim is Algorithm 3.3 which implements this idea, at least 
in the homogeneous case. The homogeneity assumption cannot be skirted, unfortunately. 
It is true that an ideal is binomial if and only if its homogenization is binomial [8, 
Corollary 1.4], but the homogenization is not accessible without a Grobner basis. It 
would be superb for our purposes if homogenizing the generators of a binomial ideal 
would always yield a binomial ideal. Unfortunately this is not the case as Example 4.1 
shows. 


3. The homogeneous case 

If a given ideal / is homogeneous, the graded vector space structure of the quotient 
k[xi,..., Xn\/I allows one to check binomiality degree by degree. For this we need some 
basic facts about quotients modulo binomials (see [8, Section 1] for details). Any set of 
binomials B in k[xi,..., x„] induces an equivalence relation on the set of monomials in 
k[xi,..., x„] under which mi ~ m 2 if and only if mi — Am 2 G {B) for some non-zero A G 
k. As a k-vector space the quotient ring k[xi,..., x„]/ {B) is spanned by the equivalence 
classes of monomials and those are all linearly independent [8, Proposition 1.11]. If the 
binomials are homogeneous, then the situation is particularly nice. For example, the 
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equivalence classes are finite and elements of the quotient have well-dehned degrees. 
The notions of monomial, binomial, and polynomial are extended to the quotient ring. 
For example, a binomial is a polynomial that uses at most two equivalence classes of 
monomials. The unihed mathematical framework to treat quotients modulo binomials 
are monoid algebras, but we refrain from introducing this notion here. 

As a consequence of the discussion above, a polynomial system /i = • • • = /^ = 0 
can be considered modulo binomials, and the coefficient matrix of the quotient system 
is well-defined. It arises from the coefficient matrix of the original system by summing 
columns for monomials in the same equivalence class. 

Example 3.1. In ]k[a;, y\, let b = x^—y'^. Among monomials of total degree three, and 
as well as x'^y and y^ become equal in ]k[x, y]/{h). Thus the degree three part in the 
quotient is two-dimensional with one basis vector per equivalence class. Consequently, 
the trinomial f = x^ + xy"^ + y^ maps to a binomial with coefficient matrix [2,1]. This 
matrix arises from the matrix [1,1,1,0] by summing the columns corresponding to x^ 
and as well as those for x‘^y and y^. 

The reduction modulo lower degree binomials in Example 3.1 can be done in general. 

Lemma 3.2. Let fi,...,fs G k[xi,...,Xn] be homogeneous polynomials of degree d, 
and B C k[xi,... ,x„] a set of homogeneous binomials of degree at most d. Then in 
the guotient ring k[xi,..., x„]/ {B) the ideal {fi,..., fs)/(B) is binomial if and only if 
the coefficient matrix of the images of fi,, fs in k[xi,..., Xn]/{B) has a partitioning 
kernel basis. 

Proof. The graded version of Nakayama’s lemma (see Corollary 4.8b together with Ex¬ 
ercise 4.6 in [9]) implies that the ideal (/i,...,/«)/ (B) C k[xi..., x„]/ (B) has a well- 
defined number of minimal generators in each degree. Therefore any minimal generating 
set consists only of degree d polynomials and Proposition 2.5 applied to the finite- 
dimensional vector space of degree d polynomials correctly decides binomiality. □ 

Lemma 3.2 is the basis for the following binomial detection algorithm. 

Algorithm 3.3. 

Input: Homogeneous polynomials /i,..., /* G k[xi,..., x^j. 

Output: Yes and a binomial generating set of {fi,..., fs) if one exists, No otherwise. 

1) Let 

• B :=(/), 

• R := k[xi,.. .,Xn], 

. F:={/i,...,/J. 

2) While F is not empty, 

a) Let Fmin be the set of elements of minimal degree in F. 

b) Redefine F := F \ F^in- 

c) Compute the reduced row echelon form A of the coefficient matrix of F^m- 
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d) If A has a row with three or more non-zero entries, output No and stop. 

e) Find a set B' of binomials in k[a;i,..., a;„] whose images in R generate 
{Fmin) and redefine B := B U B'. 

f) Redefine R := k[xi,... ,Xn]/(B). 

g) Redefine F as its image in R. 

3) Output Yes and B. 

Proof of correctness and termination. Termination is obvious. In fact, the maximum 
number of iterations in the while loop equals the number of distinct total degrees among 
/i,..., fs- Step 2.d relies on Proposition 2.2. In step 2.e, binomials that generate (Fmin) 
in R can be read off the reduced row echelon form via Proposition 2.2. Then any 
preimages in k[a:i,... ,x„] of those binomials suffice for B'. Lemma 3.2 shows that the 
while loop either exhausts F if (/i ,..., fg) is binomial, or stops when this is not the 
case. □ 

Remark 3.4. In the homogeneous case there is a natural choice of hnite-dimensional 
vector spaces to work in: polynomials of a hxed degree. In each iteration of the while 
loop in Algorithm 3.3, the rows of A span the vector space of polynomials of degree d in 
the ideal (modulo the binomials in {B)). In the general inhomogeneous situation extra 
work is needed to construct a suitable hnite-dimensional vector space. In particular, one 
needs to select from the inhnite list of binomials in the ideal not too many, but enough 
to reduce all given polynomials to binomials whenever this is possible. An interesting 
problem for the future is to adapt one of the selection strategies from the F4 algorithm 
for Grobner bases [10] for this task. 

Remark 3.5. Coefficient matrices of polynomial systems are typically very sparse. An 
efficient implementation of Algorithm 3.3 has to take this into account. 

Remark 3.6. Algorithm 3.3 could also be written completely in the polynomial ring 
without any quotients. Then in each new degree, one would have to consider the 
coefficient matrix of Fm^n together with all binomials of degree d in the ideal (B). This 
list grows very quickly and so does the list of monomials appearing in these binomials. 
Thus it is not only more elegant to work with the quotient, but also more efficient. 

To implement Algorithm 3.3 completely without Grobner bases some rehnements are 
necessary. Simply using R = k[xi,... ,Xn\/B in Macaulay2 will make it compute a 
Grobner basis of B to effectively work with the quotient. For our purposes, however, 
this is not necessary. 

Proposition 3.7. Algorithm 3.3 can he implemented without Grobner bases. 

Proof. The critical step is 2.g, when the algorithm reduces F modulo the binomials 
already found. For the following step 2.c elements of Fmin need to be written in terms 
of a basis of the hnite-dimensional vector space Fdeg(Fn,in) of degree deg(Fmin) monomials 
modulo the binomials in B. The equivalence relation introduced in the beginning of this 
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section can also be thonght of as a graph on monomials, and thus these reductions can 
be carried out with graph enumeration algorithms like breadth hrst search. Restricting 
to monomials of degree deg(Fmin), the connected components are a vector space basis 
of Rdeg(Fn,in) thus be used to gather coefficients in step 2.g. □ 

Remark 3.8. The feasibility of graph-theoretic computations in cases where Grobner 
bases cannot be computed has been demonstrated in [15]. Example 4.9 there contains 
a binomial ideal whose Grobner basis cannot be computed, but whose non-radicality 
was proved using a graph-theoretic computation. This yielded a negative answer to the 
question of radicality of conditional independence ideals in algebraic statistics. 

Remark 3.9. Using Grobner bases one represents each connected component of the 
graph defined by (B) by its least monomial with respect to the term order. Our philos¬ 
ophy is that this is not necessary: one should work with the connected components per 
se. Why bother with picking and Ending a specific representative in each component 
if any representative works? In an implementation one could choose a data structure 
that for each monomial stores an index of the connected component it belongs to. 

Remark 3.10. It is trivial to generate classes of examples where Grobner bases meth¬ 
ods fail, but Algorithm 3.3 is quick. For example, take any set of binomials whose 
Grobner basis cannot be computed and add any polynomial in the ideal. Algorithm 3.3 
immediately goes to work on reducing the polynomial modulo the binomials, while any 
implementation of Grobner bases embarks into its hopeless task. 

Remark 3.11. Remark 3.10 highlights the Grobner-free spirit of our method. The 
Grobner basis of an ideal contains much more information than binomiality. One should 
avoid expensive computation to decide this simple question. 

4. Heuristics for the inhomogeneous case 

The ideals one encounters in chemical reaction network theory are often not homoge¬ 
neous, so that the results from Section 3 do not apply. The first idea one may have for 
the inhomogeneous case is to work with some (partial) homogenization. Grobner bases 
are quite robust in relation to homogenization. For example, to compute a Grobner 
basis of a non-homogeneous ideal it suffices to homogenize the generators, compute a 
Grobner basis of this homogeneous ideal, and then dehomogenize. Although the inter¬ 
mediate homogeneous ideal is generally not equal to the homogenization of the original 
ideal, the dehomogenized Grobner basis is a Grobner basis of the dehomogenized ideal 
[2, Exercise 1.7.8]. 

Unfortunately the notion of binomiality does not lend itself to that kind of tricks. 
Geometrically, homogenizing (all polynomials in) an ideal yields the projective closure 
and dehomogenizing restricts to one affine piece. Homogenizing only the generators 
creates extra components at infinity and these components need not be binomial. Even 
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if they are binomial, the intersection need not be binomial (see also [14, Problem 17.1]). 
This is the case in the following example. 

Example 4.1. The ideal {ah — x,ah — y,x + y + l) C k[a, 6 , x, y] is binomial as it equals 
{2i/ +1, 2a; +1, 2ab+l). Homogenizing the generators, however, yields the non-binomial 
ideal {ah — xz, ah — yz, x + y + z). 

We now present some alternatives that do not give complete answers but are quick to 
check. They can be applied before resorting to an expensive Grobner basis computation. 

The quickest (but least likely to be successful) approach is to try linear algebraic 
manipulations of the given polynomials. Equivalently one applies row operations to the 
coefficient matrix, for instance, computing the reduced row echelon form. If it has a 
partitioning kernel basis, then the ideal is binomial and all non-binomial generators are 
k-linear combinations of the binomials. While it may seem very much to ask for this, it 
does happen for the family of networks in [18, Section 4]. 

If just linear algebra is not successful, one can homogenize the generators and run 
Algorithm 3.3. If the resulting homogeneous ideal comes out binomial, then the orig¬ 
inal ideal was binomial by the following simple fact, proven for instance in [ 2 , Corol¬ 
lary A.4.16]. 

Proposition 4.2. Let I C k[a;i,... ,Xn] be an ideal and I' C k[xo,xi,... ,a;n] the ho¬ 
mogeneous ideal generated by the homogenizations of the generators of I (using variable 
Xq). Then I is generated by the dehomogenization of any generating set of I'. 

We now illustrate a phenomenon leading to failure of the above heuristics. 

Example 4.3. Consider the network from [18, Example 3.15]. The steady states are 
non-negative real zeros of the following polynomials. 

fl = -kl 2 Xl + k 2 lX 2 - kui2XiX7 + {ki 2 ll + ki2l3)xg, 
f2 = ki2Xi - k2lX2 - k23X2 + ^323^3 + kgrXQ, 

/s = k23X2 — k32X3 — ^ 34 X 3 — fcggXsXy -|- kgigXg -f /CggXg, 

/4 = ^ 342^3 ~ k^QX 4 X 5 -\- kg^XQ, 

fb = —k^QX^X^ -f kg^XQ -1- kgioXg -+■ /Cl2132^95 

/e = k^QX^x^ — (^65 + kQ7)xQ, 

fl = k^jXQ — kiii2XiX7 — kggX3X7 kggXg -f A:i2iia;g, 

/s = kggX 3 X 7 — (fcgio -|- kgg)xg, 

/g = ^11123^1^7 ~ (^1211 + ^ 1213 ) 3 ^ 9 - 

The binomials /g, /g, and fg can be used to eliminate one of every pair {xq,X 4 X 5 ), 
{x 8 ,X 3 X 7 ), and (xg,XiXy). We eliminate x^x^^xg, and Xg. It can be checked that elimi¬ 
nating 0 : 10:7 instead of xg does not lead to binomials immediately (although it leads to 
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linear trinomials). 


fl = -kuXl + k2lX2, 

f2 = ki2Xi - {k21 + k 2 ‘i)x 2 + k:i 2 X 2 . + ^ 673 ^ 6 , 
/s = k2sX2 — (^32 + ^ 34 ) 0 : 3 , 

= k^/^X2, — kg^XQ, 


/5 



^1112^1213^13^7 , ksgkgioXsXr 

-keiXQ + -—-+ -——, 

^1211 + ^1213 ^910 + kgs 

, _ (^1112^1213(^910 + kgs)xi + (fci2ll + ^1213)^89^910313)3^7 

(^1211 + ^1213) (^910 + kgs) 


Using the linear relations f[ and /g the remaining system is recognized to consist of 
only two independent binomials: 


/2 

r: 

n 

/" 


—ksAXs + k^jXQ 
ksAXs — kQjXQ, 


—ke^XQ + 
kejXg + 


^ 1112 ^ 1213^21 
^ 12(^1211 + ^1213, 
^ 1112 ^ 1213^21 


+ 


^23^89^910 


(^32 + fc34)(^910 + kgs) 
^23^89^910 


^ 12(^1211 + ^ 1213 ) {ks2 + ^34)(fc910 + kgs) 


X 2 X 7 , 


X2X7. 


This analysis shows that the steady state ideal under consideration equals the binomial 
ideal {f[, f2 , fs, , fe, fs, fd) ■ The Grobner basis computation in [18, Example 3.15] 
also yields the result, but it is arguably less instructive. Note also that naive homog¬ 
enization does not yield binomiality. The element /2 is linear. After homogenization. 
Algorithm 3.3 would pick only this element as TAin and stop since it is not a binomial. 


The effect in Example 4.3 motivates our hnal method: term replacements using known 
binomials. We expect this to be very useful in applications from system biology for the 
following reasons. 

• It often happens that non-binomial generators are linear combinations of bino¬ 
mials as in Example 4.3 where fi = f[ + fg- 

• Steady state ideals of networks with enzyme-substrate complexes always have 
some binomial generators. These complexes are produced by only one reaction 
and thus their rate of change is binomial. 

• In MAPK networks, which describe certain types of cellular signaling, one often 
hnds binomials of the form kXaXh — k'xc- 

• Frequently binomials in steady state ideals are linear. Equivalently some of 
the concentrations are equal up to a scaling (which may depend on kinetic 
parameters). This happens for all examples in [19]. 
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We now illustrate term replacements in a larger example which comes from the net¬ 
work for ERK activation embedded in two negative feedback loops (see [5, Section 5] 
for pointers to the relevant biology). 

Example 4.4. Consider the following steady state ideal generated by 29 polynomials. 

/i = -kiXiX 2 + k2X3 + keXe, f 2 = -kiXiX 2 + fexa + * 3 X 3 , fs = kiXiX 2 - ^ 2 X 3 - ^ 3 X 3 , 

/4 = fciiXio + fcl 2 Xio + k3sX25 + ki2X27 + kzX 3 - ^ 37 X 13 X 4 , 

— ^ 4 X 4 X 5 + fcsxe — ^ 7 X 4 X 7 + /cgxs + fcgxs — ^ 10 X 4 X 9 
/s = A:i4Xi 2 + fci5Xi2 + fci7a;i3 + ki^xiz + ^ 35 X 24 + ^ 35 X 24 + ^ 41 X 27 + ^ 42 X 27 - ^ 13 X 11 X 5 - ^ 34 X 16 X 5 
^^ 40 X 29 X 5 A^ 4 X 4 X 5 “t” k^XQ -\- k^XQ /tigXaXg, 

fe = ^4x4x5 - fcsxe - fcexe, fr = kisxis — ^7x4x7 + ksxs, /s = kjx^x^ - ksxs - kgxs, 
h = ^ 112^10 + fci5a;i2 + fci7a;i3 + k^x^ — ^ 10 X 4 X 9 — kiQX^xg, fio = —knxio — kigxio + ^ 10 X 4 X 9 , 
fll = kl2Xl0 + fcl4Xl2 — fcl9XllXl4 + feoXls + felXls — fegXnXie + k23Xl7 + ^24X17 — /ciaXuXs, 

/l 2 = —^142^12 — fcl52;i2 + Al 32 ;ilX 5 , /13 = —^ 17 X 13 — fclsXia + kiQX^Xg, 

flA = —fcl92;ilXi4 + fc202;i5 + ^302^21 + k3QX24, fl5 = fcl92;ilXi4 — ^ 29 X 15 — ^ 21 X 15 , 

fl6 = k2lXi3 — k22XllXie + k23Xn — k2SXieXiQ + k27X20 + ^ 92^21 + ^ 332^23 + ^ 352^24 ~ fc342^162;5j 

fl7 = ^ 222 : 11 X 16 — ^232^17 ~ ^ 42 : 17 , 

/18 = ^242^17 ~ k 25 XisXig + k 2 &X 20 — ^312^182^22 + ^322^23 + ^382125 + ^392^25 

— A:43Xi8X28 + k44,X2Q + ^453^29 ~ k37XisX4, 

fig = —k4eXig — kgsXigXig — ^ 25 X 13 X 19 + kgeXgo + ^72120 + kggXgi + k3gX21 + fc452l29j 

/20 = fc252;i321l9 ~ kgeXgg — ^72120, / 2 I = A:232;i621l9 — kggXgi — A:3oX2l, 

/22 = —k3lXi8X22 + ^22123 + ^332123, /23 = fc3l21l32;22 — ^322^23 “ ^332123, 

/24 = —k35X24 — k3eX24 + ^34X15X5, fgb = —^382125 ” ^392^25 + k 37 Xi 8 X 4 , 

/26 = ^392125 + ^4ia;27 “ ^402l262l5) /27 = —^4ia;27 ” ^422127 + fc402;26a;5 

/28 = fc462;i9 “ A:43 Xi 3X28 + A:44X29, /29 = ^432;i82l28 ~ ^442129 “ k43X2g- 

After some obvious factorization, the following elements are binomials: / 2 , fe, fg, fio, 
fi 2 , fi 3 , fi 5 , fn, / 20 , / 21 , / 22 , / 23 , / 24 , / 25 , / 27 , / 29 - The system has seven conservation 
relations, which can be found by linear algebra. According to our strategy to use 
binomials to simplify the system, we eliminate, if possible, non-binomials using the 
conservation relations. This is not always possible, as some of the conservation relations 
stem from duplicate equations like /2 = —fg. We eliminate /s, Ja, /s, fg, /g, fig, and fig. 
The remaining non-binomial part consists of /i, /y, fn, fu, /le, / 26 , and f 2 g. Dividing 
by reaction constants, each of the binomials is of the form Xi = KxjXi for some rational 
expression K involving only reaction constants. Using these in the non-binomials yields 


f'i = 
/n = 


kik 3 XiX 2 kAkgXiXz 

kg -1-^3 /cs + fee ’ 

^ 13 ^ 152 : 11 X 5 kigki 2 X 4 Xg 
ki4 + ki3 kii + ki2 


f7 = 


k7kgX4X7 ^ kiekisX^Xg 
kg + kg ki7 + kis 
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fU = 


fl6 — 


/26 ~ 


fcl9feia^lia^l4 ^28^302^162^19 

/C20 + ^21 feg + ^30 

^19fel2;il2;i4 fc22fc242;il2:i6 

^20 + fel ^23 + ^4 

, k3ik33XisX22 , 

H-;--A;34a;i62;5 + 

K 32 + ^33 

^37^392;i8a;4 ki3ki2X2QX3 
kss + ^39 ^41 + ^42 


^ ^34^362'162'5 
^35 + ^36 

k2sk29XieXig 
k28Xl6Xig H- - --- 

K 29 + kso 

kzik^^XiQX^ 

^35 + ^36 


^25^272^182^19 
^26 + ^27 


/28 — ^462^19 


kii3ki3Xi3X28 
A^44 + /C45 


In particular, we find five new binomials /(, /y, f[^, /^g, and /^g. Adding /{4 to /{g 
yields the trinomial 


J./, _ ^22^24Tlia;ig ^ fc25^272^182;i9 ^31^332^182^22 

^23 + ^24 ^26 + ^27 ^32 + ^33 

Consequently, the original system is equivalent to a system consisting of 27 binomials 
and two trinomials of a relatively simple shape. For comparison we computed the 
Grobner basis in Macaulay2 with rational functions in the reaction rates as coefficients. 
Although the computation hnished in just 18 minutes, the result is practically unusable. 
The Grobner basis consists of 169 elements each of it with huge rational functions as 
coefficients. The structure that we observed above is completely lost. 


The lesson learned from Example 4.4 is that term replacements using binomials are 
useful in solving a polynomial system, even if the end result is not binomial. Especially 
in the non-homogeneous case where the notion of minimal generators is absent, compu¬ 
tations with the Binomials package [16] in Macaulay2 [11] can probably only assist, 
but not automatically do useful reductions. For example, a natural general choice would 
be to replace higher degree monomials by lower degree ones, but this would not directly 
reveal binomiality in Example 4.3. 

Finally, we summarize a possible strategy to deal with inhomogeneous ideals. Exam¬ 
ple 4.1, for instance, is solved already by item 1, but also by item 3. 


Recipe 4.5. 

1) Try linear algebra and Proposition 2.5. 

2) Homogenize the given ideal and run Algorithm 3.3. If the algorithm returns bi¬ 
nomials, then by Proposition 4-3 the original (dehomogenized) ideal is binomial. 
The homogenization should be earned out after linear algebra reductions to pos¬ 
sibly detect homogeneity already at an earlier stage (compare Example 2.8). 

3) If Algorithm 3.3 returns a negative answer, dehomogenize and use known bino¬ 
mials for term replacements (as in Example 4-4) ■ Potentially homogenize again 
with an enlarged generating set. 

4) Compute a reduced Grobner basis. 
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